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C/^ ; ABSTRACT 

Some merging galaxy clusters host diffuse extended radio emission, so-called radio halos and relics, unrelated to individual galaxies. 
The origin of these halos and relics is still debated, although there is compelling evidence now that they are related to cluster merger 
events. Here we present detailed Westerbork Synthesis Radio Telescope (WSRT) and Giant Metrewave Radio Telescope (GMRT) 
radio observations between 147 MHz and 4.9 GHz of a new radio- selected galaxy cluster IRXS J0603. 3+4214, for which we find 

Ua redshift of 0.225. The cluster is detected as an extended X-ray source in the ROSAT All Sky Survey with an X-ray luminosity of 
Lx, 0.1-2.4 kev ~ 1 X 10^^ erg s~^ The cluster hosts a large bright 1.9 Mpc radio relic, an elongated ~ 2 Mpc radio halo, and two fainter 
, smaller radio relics. The large radio relic has a peculiar linear morphology. For this relic we observe a clear spectral index gradient 

O I from the front of the relic towards the back, in the direction towards the cluster center. Parts of this relic are highly polarized with a 

I ■ polarization fraction of up to 60%. We performed Rotation Measure (RM) Synthesis between 1.2 and 1.7 GHz. The results suggest 

O ' that for the west part of the large relic some of the Faraday rotation is caused by ICM and not only due to galactic foregrounds. 

I We also carried out a detailed spectral analysis of this radio relic and created radio color-color diagrams. We find (i) an injection 

C/^ . spectral index of -0.6 to -0.7, (ii) steepening spectral index and increasing spectral curvature in the post-shock region, and (iii) an 

■ overall power-law spectrum between 74 MHz and 4.9 GHz with a = -1.10 + 0.02. Mixing of emission in the beam from regions 
with diff'erent spectral ages is probably the dominant factor that determines the shape of the radio spectra. Changes in the magnetic 
field, total electron content, or adiabatic gains/losses do not play a major role. A model in which particles are (re) accelerated in a 
first order Fermi process at the front of the relic provides the best match to the observed spectra. We speculate that in the post-shock 
region particles are re- accelerated by merger induced turbulence to form the radio halo as the relic and halo are connected. The 
orientation of the bright relic and halo indicate a north- south merger event, but the peculiar linear shape and the presence of another 
relic, perpendicular to the bright relic, suggest a more complex merger event. Deep X-ray observations will be needed to determine 
the merger scenario. 
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1, Introduction ically unpolarized and have an extent of about a Mpc. 

^ , , There are two main classes of models proposed to ex- 

Radio halos and rehcs are diffuse sources found in galaxy clus- ^j^^ ^^^^ j^^l^^. relativistic electrons in 

. to^t hat show s igns of a n ongoing m erger event (see Feretti ^j^^ j^M are re-accelerated in-situ through interaction with 

^ • mi Ferraneta l. 2008; Bru£geneUL2011, for an overview), turbulence gen erate d in the I CM by cluster-cluster mergers 

- ■ Ihese sources reveal tne presence oLremstic particles and terunetti et al.ll200H lPetrosianll2001h . or (2) the energetic elec- 

cluster-wide magnetic fields (e.g., MimW within the intr- ^^^^^ secondary products of proton-proton collisions (e.g., 

acluster medium (ICM). In the hierarchical model of structure fDennison 1980; Blasi & Colafrancesco 1999; Dolag & EnBliiJ 

formation galaxy cluster form through a sequence of mergers 200O; MiniatietaU 120011; iKeshet & Loebl I2010t lEnBlinetall 

with snia Her substructures^ These merger events create shocks ^^^^-^^^ ^^^^j^ ^ combination of th e two 



X: 



and turbulence in the ICM and could amplify magnetic fields, n^echanisms hav e 'been considered (iBrunetti & Blasil 120051 

Several models have been put forward which link the pres- Ippiag et al.ll200alBrunetti & Lazarianll20 1 1) 
ence of these diffuse radio sources to clust e r merger events 

re,g.. lEnBlin et al." 7998": 'Brunetti et al." •2001: "Petrosian '200il ^^^^o rg//c^ are usually divided up into three different 

&sano & Brunetti 2005; Keshet 2010; EnBlin et al. 201 1). classes (see| Kempneret alj | 2004 for an overview). (1) Radio 

Radio halos are smooth extended radio sources that roughly <^'^c/!f are large extended sources mostly found in the outskirts 

follow the X-ray emission from the ICM. They are typ- of galaxy clusters They are proposed to trace shock waves 

( En Blin et all 1 19981; iMiniatiet all I200lh . in which particles are 

* Einstein Fellow, e-mail: rvanweeren@cfa.harvard.edu accelerated by the diffusive shock acceleration mechanism 
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(PSA: e.g.. rKmnskii|fT977l: lAxfor^ et aDll97Zl; IMillQTSallbl: 
'Blandford & Qstriked fT978l: iDrurv "1983; "Blandford & Eichler| 



2. Observations & Data Reduction 



19911: iMalkov & O'C Drury 2001) in a 2.1. GMRT observations 



,1987; Jones & Ellison 
first-order Fermi process. However, the efficiency with which 
collisionless shocks can accelerate particles is unknown and 
may not be enough to produce the observed radio brightness of 
relics. A closely linked scenario is that of shock re- acceleration 
of pre-accelerated electrons in the ICM, which is a moreeffi- 
cient mechanism for we ak shocks (e.g., Markevitch et al.ll2005l: 
iGiacintucci etaP 120081: iKang & Rvui l2011: Kang etal.1 l2012h . 
(2) AGN relics trace old radio plasma from previous episodes 
of AGN activity. (3) Radio phoenices are AGN relics that 
have been compressed adiaba tically by merger shock waves 
boosting the radio emi ssion (lEnBlin & Gopal-Krishnal 120011: 
lEnBlin & Briiggenl l2002h . The radio spectra of these sources 
should be steep {a < -1.5, oc y^, where a is the spectral 
index) and curved due to synchrotron and Inverse Compton (IC) 
losses. 

A diff'erent model for radio rel ics (radio gischt) and halos 
has been proposed by lKeshej (12010 ). It is based on a secondary 
cosmic ray electron model, where the time evolution of the mag- 
netic fields and cosmic ray distribution are taken into account 
to explain both halos and gischt. For the outer edges of relics 
it is predicted that the spectral index should be a ^ - 1 . In ad- 
dition, the magnetic fields {B) are expected to be strong, with 
B > BcMB > 3 yuGauss, and ^cmb the equivalent magnetic field 
strength of the cosmic microwave background. 

Although there is now substantial evidence that radio 
relics and ha l os ar e related to galaxy cluster mergers (e.g., 
ICassano et al.l l201Qh . the detailed physics are still not under- 
stood. Radio spectra are a crucial way to separate between 
the physical models for the origin of radio relics and halos. 
Obtaining high-quality spectra is however difficult as these 
sources are extended and have a low surface brightness. In par- 
ticular, spectral curvature provides important information about 
the underlying physical processes. However, curvature in the 
synchrotron spectrum is only observable if the observed fre- 
quency range is wide enough. In addition, obtaining images at 
matched resolutions over a wide-frequency range is often not 
possible. 

We carried out an observing c ampaign to searc h for new 
diff'use radio sources in clusters (van Weeren et al. ^ 2011c,b). 
By inspecting the W ENSS (Rengelink et al.lll997h and NVSS 
(ICondon et al 1 1 19981) surveys we came across the complex ra- 
dio source B3 0559-F422B. The radio source roughly coincided 
with an extended X-ray source in the ROSAT All Sky Survey. In 
addition, an overdensity of galaxies following the X-ray emis- 
sion, was visible in 2MASS images. The diff'use radio source, 
extended X-ray emission and galaxy overdensity strongly sug- 
gested the presence of a previously unidentified galaxy cluster 
located at moderate redshift. In this paper we present detailed ra- 
dio observations of this new galaxy cluster IRXS J0603.3-h4214 
and investigate the spectral and polarimetric properties of the 
diff'use emission. 

The layout of this paper is as follows. In Sect. [2] we give an 
overview of the observations and the data reduction. The WSRT 
and GMRT images are presented in Sect. [3l The radio spectra 
and polarization data are analyzed in Sects |4]to[6l We end with a 
discussion and conclusions in Sects. |7] and [H 

Throughout this paper we assume a ACDM cosmology with 
Hq = 71 km s"^ Mpc"\ = 0.3, and Qa = 0.7. All images are 
in the J2000 coordinate system. 



GMRT obs ervations \yere ta ken using the GMRT software back- 
end (GSB; iRov et al.ll2QTQh . An overview of the observations is 
given in Table [T] 

The data were reduced with the NRAO Astronomical Image 
Processing System (AIPS) package. The data were visually in- 
spected for the presence of radio frequency interference (RFI) 
which was subsequently removed. For the 147 and 241 MHz 
data, RFI was fitted and subtracted using the techn ique de - 
scribed bv lAthreval (l2009l) . implemented in Obit (.Cotto 



iQue de - 
inl l2008h . 



Standard bandpass and gain calibration were carried out, fol- 
lowed by several rounds of phase self-calibration and two final 
rounds of amplitude and phase self-ca libration. The flu xes for 
the calibrator s ources were set by the IPerley & Tayloj (Il999b 
extension to the'Baars et al. scale. Images were made us- 

ing "briggs" weighting (Bri ggsll 19951) . see Table[T] were cleaned 
down to 2 times the rms noise level (2crj:ms) within the clean 
boxes and corrected for the primary be am attenuatiorQ. For mor e 
details about the data reduction see ivan Weeren et al.l (l2011dh . 
The 147 MHz data were further calibrated for ionospheric phase 
distortions, as these can be come quite severe a t this frequency, 
with the SPAM package (llntema et al.l '2009). At 325 MHz 
we removed se veral sources using the "peeling"-method (e.g., 
lNoordamll2Q04h . We assu me a 5% uncertainty in the calibration 
of the absolute flux-scale (iChandra et al.ll2004l) . 



2.2. WSRT observations 

WSRT observations of IRXS J0603.3-h4214 were taken in the 
L-band, 13cm and 6cm bands, see Table [2] for details. The data 
were calibrated with the CAS^aEI package. We first removed time 
ranges aff'ected by shadowing and RFI. The data were then cal- 
ibrated for the bandpass response and subsequent gain solutions 
were transferred to the target source. Channel dependent leakage 
terms were found using an unpolarized calibrator source and the 
polarization angles were determined from 3C2 86. The fluxes for 
the calibrator s ources were set by the Perley & Taylod ([l999) ex- 
tension to the lBaars et al.l (Il977h scale. The data were exported 
into AIPS for several rounds of phase self-calibration, followed 
by two rounds of amplitu de and phase self-calibration, see also 
van Weeren et al. (l2QlTdh . Images were cleaned with manually 
placed clean boxes and corrected for the primary beam attenua- 
tion. A deep image was created by combining the images from 
individual IF's from the 18, 21, and 25 cm bands. The images 
were convolved to the same resolution and combined with a 
spectral index scaling of -1. The resolution of this combined 
image is 29'' x 19''. We assume a 5% uncertainty in the cali- 
bration of the absolute flux-scale. 

We used the WSRT L-band observations to perform RM- 
Synthesis (see Sect. 16.2b . For this, we created full polarization 
maps of every two neighboring frequency channels (i.e., the im- 
age bandwidth is 62.5 kHz). We inspected all these images and 
removed maps that had high noise levels, were afl'ected by RFI 
or had other artifacts. In the end, 280 channels maps between 
1 170 and 1786 MHz were retained for RM-synthesis. 



^ http://gmrt.ncra.tifr.res.in/gmrt_hpage/Users/doc/manual/ 
User sManual/node27. html 
^ http://casa.nrao.edu/ 
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Table 1. GMRT observations 





147 MHz 


241 MHz (dual) 


325 MHz 


610 MHz (dual), 610 MHz 


1280 MHz 


Observation date 


Apr 29, 2010 


Nov 20, 2009 


Apr 25, 2010 


Nov 20, 2009 & 30 Apr, 2010 


Nov 12, 2009 


Usable bandwidth 


ij iVlrlZ 


D iVlrlZ 


JU iVlrlZ 


JU iVlrlZ 


jU MrlZ 


Channel width 


31.25 kHz 


62.5 kHz 


62.5 kHz 


62.5 kHz 


62.5 kHz 


Polarization 


RR+LL 


LL 


RR+LL 


RR, RR+LL 


RR+LL 


Integration time 


4 sec 


8 sec 


8 sec 


8 sec 


16 sec 


Total on- source time 


5.5 hr 


3.5 hr 


6.5 hr 


3.5 + 6.0 hr 


4.0 hr 


Beam size 


26" X 22" 


16" X 12" 


12" X 8.7" 


5.1" X 4.1" 


5.8" X 2.0" 


Rms noise ((Trms) 


0.92 mJy beam"^ 


300 yL/Jy beam~^ 


79 yL/Jy beam~^ 


26 yL/Jy beam~^ 


32 yL/Jy beam~^ 


Briggs weighting 


0.5 


0.5 


0.0 


0.0 


0.5 



Table 2. WSRT observations 



25 cm, 1221MHz 21 cm, 1382 MHz 18 cm, 1714 MHz 13 cm, 2272 MHz 6 cm, 4.9 GHz 

Bandwidth 8 x 20 MHz 8x20 MHz 8x20 MHz 8 x 20 MHz 8x20 MHz 

Number of channels per IF 64 64 64 64 64 

Polarization XX, YY, XY, XY XX, YY, XY, XY XX, YY, XY, XY RR, LL, RL, LR XX, YY, XY, XY 

Observation dates Aug 28, 2010 Sep 10, 2010 Sep 3, 2010 Sep 9, 2010 Sep 2 &6, 2010 

Total on- source time 12hr 12hr 12 hr 12 hr 24hr 

Beam size 29.0" x 19.0" 27.1" x 16.4" 23.6" x 15.6" 16.1" x 9.8" 7.0" x 4.7" 

Rms noise (cTj-ms) 39 fiJy beam~^ 27 yt/Jy beam~^ 25 yt/Jy beam~^ 40 yt/Jy beam~^ 41 yt/Jy beam~^ 

Briggs weighting 0.0 0.5 1.0 0.5 0.5 



2.3. WHT spectroscopy & imaging 

Optical images of IRXS J0603.3+4214 were taken with the 
PFIP camera on the 4.2m William Herschel Telescope (WHT) 
between 15 and 19 April, 2009 with V, R and I filters. The see- 
ing varied between 1.0'' and 1.5''and the total integration time 
was 1500 s per filter. The data were reduced with IRAF (iTodvl 
Il986lll993l) and the mscred package ( Valdes 1998). Ima Res were 
flat-fielded and corrected for the bias ofl'set. R and I band images 
were also fringe corrected. The individual exposures were aver- 
aged, with pixels being rejected above 3.0crrms to remove cos- 
mic rays and other artifacts. The optical images were strongly 
afl'ected by the bright star BD+42 1474 (Vmag^ 8.62), see ¥ig.\f\ 
To determine the redshift of the cluster, WHT ISIS spectra 
of galaxies located on the optical images were taken on 10 and 
11 February, 2011. For the blue arm we used the R300B grat- 
ing and for the red arm the R316R grating. The slit- width was 
1.5''. Flat-fielding, bias correction, and wavelength calibration 
were performed in IDL0. The total exposure time per galaxy was 
1500 s in both the blue and red arms. In the end, the spectra for 
the blue and red arm were merged into single spectrum for each 
galaxy. 



3. Results: redshift, X-rays, and radio continuum 
maps 

3.1. Redshift of 1 RXS J0603.3+4214 

The brightest cluster members are visible in the 2MASS survey 
images ( Skrutskie et al. 2006) and some of them are also listed 
in the 2MASS Extended sources catalog (Skruts kie et al.ll2003!i) . 
The cluster is not listed in NED or SIMBAD and probably re- 
mained unidentified because it is located relatively close to the 
galactic plane at a galactic latitude b = 9.4° and nearby the bright 
star BD+42 1474. 



^ http://www.ittvis.com 



Table 3. Redshifts 



Galaxy 


2MASS K mag z 


2MASX J06031667+4214416 


13.204 


0.227 


2MASX J06030757+4216215 


13.610 


0.222 


2MASS J06032605+4214050 


15.075 


0.228 


2MASS J06030644+42 15241 


15.201 


0.227 


2MASX J06032432+4209306 


13.858 


0.220 



No published redshifts are available for any of the galax- 
ies in the cluster. Galaxies for which WHT ISIS spectra were 
obtained are listed in Table [3] and the spectra are shown in 
Fig. [21 The spectra are typical for passive red elliptical galax- 
ies mostly found in clusters. We find that the five galaxies are 
located at 0.220 < z < 0.228. Taking the average value we adopt 
z = 0.225 ± 0.04 for the cluster, with the uncertainty in the red- 
shift given by the standard deviation. 

3.2. X-ray emission from the \Cf\A 

IRXS J0603. 1+4214 is seen the ROSAT All-Sky Survey as 
an extended source and fisted as IRXS J0603 13. 4+42 1231, 
IRXS J060322.3+4213 05, and IRXS J0603 14.8+421439 
(IVoges etal.lll999L 120001) . see Fig. [3] (left panel). Using the red- 
shift and ROSAT count rate (0.21 PSPC cts s"^) we find an X-ray 
luminosity (Lx, 0.1-2.4 kev) of ~ 1 x 10^^ erg s"^ . The X-ray emis- 
sion is extended in the north- south direction. Additional emis- 
sion extends to the east and west of the main X-ray peak. The 
high X-ray luminosity and extended emission are consistent with 
a massive cluster undergoing a major merger event. 

3.3. Radio continuum maps 

The WSRT L-band map is displayed in Fig.O We labeled some 
of the most prominent sources in this image. The GMRT 147, 
241, 325, 610, 1280 MHz, and WSRT 21 and 13 cm images are 
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Fig. 1. WHT V, R, I color image. Yellow contours are from the GMRT 610 MHz image (Fig. [4]) and drawn at levels of [1 , 2, 4, 8, . . .] x 
0.15 mJy beam"^ Galaxies for which spectra were obtained are marked with red circles, see also Table [3] and Fig.O The spiral 
pattern (yellow contours) at the bottom right is due to residual calibration errors around the source B3 05 5 9 +422 A. 



shown in Figs. |4] to [6] (for the 4.9 GHz image see Fig.[T71). The 
properties of the diffuse source in the cluster are summarized in 
TablelH 



Source A (B3 0559+422A) is a com] 
source with a flux of 0.29 Jy at 1.4 GHz dJackson et a 



iBeaslev et all l2002t iMarecki et aP 119991: IPatnaik et al 



spectrum 
120071: 



1992). 



The radio source is associated w ith a star-like object in the WHT 
images and listed as a quasar by 'Andrei et al.l (l2009h . but no red- 
shift is reported. Source F is a compact source (but resolved in 
the 1280 and 610 MHz images) and has an optical counterpart 
in the WHT images. 

The most prominent source in the GMRT and WSRT im- 
ages is B3 0559-h422B (sourc e B). The source is also visible 
in the 74 MHz VLSS survey (ICohen et al.ll2007l) and Hsted as 



VLSS J0603.2-h4217. In the L-band image, the source has a 
largest angular size of 8.7 \ which corresponds to a physical 
size of almost 1900 kpc. The source consists of a bright west- 
ern part (Bl) and a fainter linear extension to the northeast (B2, 
B3). The radio emission brightens and fades two times along this 
extension, while the width of the source also varies. In the high- 
resolution 1280, 610, and 325 MHz images, the source displays a 
complex filamentary morphology. Some "streams" of emission 
extent from the northern part of Bl to the south. The northern 
boundary of the source is sharp, while the emission fades more 
slowly at the southern part. We classify source B as a radio relic 
because of the lack of optical counterparts, its polarization and 
spectral properties (see Sects. |4]to[6]), and the large physical elon- 
gated size and peripheral location in the cluster. 
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Fig. 3. Left: X-ray emission from ROSAT. The image from the ROSAT All Sky Survey was smoothed with a 200'' FWHM 
Gaussian and is shown in orange colors. Solid contours are from the WSRT L-band image and drawn at levels of [1, 2, 4, 8, . . .] x 
0.15 mJy beam"^. Right: WSRT L-band image (1 160-1780 MHz) with sources labelled. Contour levels are drawn at [1, 2, 4, 8, . . .] x 
80 //Jy beam"^. The beam size is 29'' x 19" and indicated in the bottom left corner of the image. 
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Table 4. Diffuse radio sources in IRXS J0603.3-h4214 



Source 


S 1382MHz 


^ 1.4GHz 


LLS^ 


a 




mJy 


10^4 


kpc 




B 


319.5 + 20.8 


60 


1870 


-1.10 + 0.02 


C 


35.9 + 2.6 


6.8 


1700-2100 


-1.15 + 0.06 


D 


5.38 + 0.31 


1.0 


215 


-1.10 + 0.05 


E 


9.48 + 0.97 


1.8 


860 


-1.0 + 0.2 



^ largest linear size 



Fig. 2. Rest- frame WHT ISIS spectra for the galaxies listed in 
Table[3l The order (from top to bottom) is the same as in Table[3] 
Various absorption features are indicated. 



Extending from B 1 to the south, there is a patch of low sur- 
face brightness emission (C: subdivided into CI and C2). A hint 
of emission is visible north of Bl in the L-band and 325 MHz 
images, possibly associated with CI. The surface brightness de- 
creases more rapidly at the south side of C (C2) than at the east 
and west sides. If we consider the faint emission north of Bl to 
be part of C then the source has a largest angular size of about 
10' (i.e., a physical extent of 2.1 Mpc). Only counting the emis- 
sion south of Bl the source extends about 8'. We classify source 
C as the giant elongated radio halo because it roughly follows the 
X-ray emission from the ICM and has different spectral index 
properties than source B. The radio power of 6.8 x 10^^ W Hz"^ 



falls on the Lx-/ ^ i.4GHz correlation for giant radio halos (e.g., 
iLiang et al.l 120001: ICassano etaDl2006h . However, at the most 
southern part of C2 the surface brightness increases slightly. This 
is best seen in the 147 and 241 MHz images (Fig. [6]). The south- 
ern end of C2 could therefore be the "counter" radio relic of Bl 
(see Sect. IT]). Although, we can also not exclude the possibility 
that it is part of the radio halo. 

Another diffuse elongated source is located to the east of CI . 
The source consist of two parts. El and E2. We classify source 
E as a radio relic because we could not find an optical counter- 
part associated with it and it is located at the eastern boundary 
of the X-ray emission. The extent of E is 4' which corresponds 
to 860 kpc. Source D is another diffuse source without any op- 
tical counterpart in the WHT images, which we also classify as 
a radio relic. It has a largest angular size of 1.0' in the GMRT 
610 MHz image and displays hints of a complex morphology 
(Fig.O, very different from a typical radio galaxy. 




Fig. 4. GMRT 610 MHz image. Contour levels are drawn at V[l. 2,4, 8, . . .] x 4(Trms- Negative -3(Trms contours are shown by the 
dotted Hnes. The beam size is S.T' x 4.r' and indicated in the bottom left corner of the image. 



4. Radio spectra 

4.1. Spectral index maps 

We created radio spectral index maps using the GMRT and 
WSRT data. We first made radio maps at roughly the same reso- 
lution applying suitable tapers in the uv-plane and using uniform 
weighting. We discarded data below 0.25 k/l to select only com- 
mon inner uvranges. The maps were then convolved to the same 
resolution. A high-resolution (1.9'' x 6.2'') spectral index map 
between 610 and 325 MHz is shown in Fig.[7l Pixels with values 
below 5(Trms wcrc blanked. A medium-resolution (20'' x 18") 
spectral index map between 2272 and 147 MHz was also made 
by fitting a second order polynomial in log (S) - log (v) space 
through the flux measurements at 2272, 1714, 1382, 1221, 610, 
325, 241, and 147 MHz, see Fig. [8] (left panel). Pixels with values 
below 1.5(Trms wcrc blanked. To map the spectral index across 
the low surface brightness halo emission in the cluster, we also 
made a low-resolution spectral index map. We convolved the 
eight maps with a 35" FWHM Gaussian (giving a resolution 
of approximately 40") and fitted a power-law spectral index to 
minimize the number of fitted parameters. Pixels with a spectral 
spectral index error > 0.7 were blanked (Fig. [HI right panel). The 
errors in the spectral index maps are computed on the basis of the 
(Trms values for the individual maps at the various frequencies 
and the reported flux calibration uncertainties in Sect. O When 
maps at more than two frequencies are used, it is assumed that 
the adopted fitting function provides a good physical descrip- 



tion of the spectral shape. Figures displaying the corresponding 
errors in the spectral maps are show in Appendix lAl 

The spectral index across relic B displays a clear north-south 
gradient, with a spectral index of about -0.6 to -0.75 on the 
north side of the relic, steepening to ~ -1.9 at the south side 
of B2 and B3 and to < -2.5 for Bl. The spectral index gra- 
dient is visible over the entire length of the relic. The high- 
resolution 610-325 MHz spectral index map reveals the same 
general trends for relic B as in the medium-resolution map. 
However, the SNR on the spectral index is somewhat lower, 
mainly because of the smaller frequency span. We note that the 
spectral index at the front of B2 and B3 has a steeper spectral in- 
dex (about -0.9 to -1.1) in the region where the surface bright- 
ness drops (i.e., at the intersection Bl-B 2 and B2-B3). 

The spectral index across Relic E varies mostly between 
-1.0 to -1.2. There are no systematic trends visible across the 
relic, except that the north part of E3 has a somewhat flatter spec- 
tral index, see Fig.[8l The spectral index for relic D steepens from 
-1.0 to -1.3 from south to north. The spectral index across the 
radio halo (C) is difficult to determine as the SNR is low, but it 
roughly varies between -1.5 and -0.8 with an average of -1.1. 
There is a hint of spectral flattening for the center of the halo, 
while to the north and south of it the spectral index is steeper on 
average. 
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IRT 150, 241, 325, 1280 MHz and WSRT 1382 and 2272 MHz images. Contour levels are drawn at [1, 2, 4, 8, . . .]x4crn 
-3(Trms contours are shown by the dotted Hnes. The beam size is shown in the bottom left corner of the images. 
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Fig. 5. GMRT 610 MHz cutout around sources D and E. 
Contours are drawn as in Fig.|4l 



4.2. Integrated radio spectra 

We determined the integrated radio spectra of the components, 
taking the same maps which were used to create the medium- 
resolution (2C X 18'0 spectral index map to minimize the ef- 
fects of different uv-coverage. In addition, we added flux mea- 
surements at 74 MHz and 4.9 GHz for reHc B from the VLSS 
survey and the WSRT observations. The radio spectrum is shown 
in Fig. [9] and is well fitted by a single power-law with a = 
-1.10 ± 0.02 between 74 and 4900 MHz. For relics D and E we 
find integrated spectral indices of -1.10+0.05 and -1.0+0.2. For 
source C we find a = -1.15 + 0.06. The integrated radio spec- 
tra do not show clear evidence for spectral breaks or turnovers, 
although the radio spectrum for E is poorly determined. 



5. Radio color-color diagrams 

We investigated the detailed spectral shape of the relic emis- 
sion. We first divided the relic into two parts, Bl and B2+B3, 
which we analyzed separately. These parts were further subdi- 
vided into regions (ofref ) within a specific spectral index range of 
0.1 units, using the spectral index map from Fig. [8] (left panel). 
The resulting spectra, including first order polynomial fits (in 
log(/) - log(y) space) through the data points are shown in 
Fig. \T3\ The fluxes were normalized by dividing by the number 
of pixels found in each region, i.e., we work in surface bright- 
ness units. Each spectra thus corresponds to a region from the 
spectral index map where a^Qf - 0.05 < a^Qf < a^Qf + 0.05, with 
a,^f = -0.65,-0.75,-0.85, etc. We limited a^^f to -1.85 and 




Frequency [MHz] 

Fig. 9. Integrated radio spectrum for radio relic B. 

-1.65 for Bl and B2+B3, respectively, to retain sufficient SNR 
at the highest frequency map (at 2272 MHz). We use these spec- 
tra as a starting point for creating the radio color-color diagrams. 

5.1. Spectral models 



In the model put forward by lEnBlin et"aD ([1998"), relics trace 
shock waves in which particles are accelerated by the DSA 
mechanism. Without projection effects and mixing of emission, 
and all properties of the shock remaining constant, the spectra 
at the different locations should simply reflect the energy losses 
of the radiating particles. Assuming an edge-on planar shock- 
wave, the time since acceleration for particles at a given loca- 
tion behind the front of the relic is simply //vd, with / the dis- 
tance from the front of the shock and Vd the shock downstream 
velocity. According to DSA-theory, the injection spectral index 
QTini is linked to the Mach number (At) of the shock by (e.g., 
iBlandford & Eichlejll987h 



1 At^ + 1 



1 



"""J "2 M2 . 



(1) 



The integrated spectral index is steeper by about 0.5 units com- 
pared to Qfinj for a simple shock model where the electr on cool- 
ing time is much shorter than the lifetime of the shock (iMiniatil 
2002). Most of the synchrotron radiation we observe comes from 
the "critical" frequency Vc 



(2) 



Directly behind the front of the shock, the spectra should have 
a power-law shape of the form I{y) = /qv^^^j , under the usual 
assumption that one starts with a power-law distribution of rela- 
tivistic electrons N(E) = NoE~\ where ^ = 1 - 2ainy 

S ynchrotron spectra can be calculated using the standard for- 
mula (|Pacholczv]<j|l970h 



/(y, t) oc 



r sin(^) r 

J47: Jo 



F(x)N(E, 0, t)dEd^l, 



(3) 



where N(E, 6, t) is the electron energy distribution, 6 the pitch 
angle between the electron velocity and magnetic field vectors, 
X = y/vc, and F{x) = Ks/3(z)dz, with K^i^iz) an irregular 
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Fig. 7. GMRT 610-325 MHz spectral index map. Contour levels from the GMRT 325 MHz image are drawn at [1,4, 16, 64, 
6(Trms and the spectral index map has a resolution of 1.9" x 6.2''. Pixels below 5(Trms are blanked. 



modified Bessel function. Synchrotron and IC (radiation) losses 
change the electron energy distribution over time at a rate 



at 



(^syn + ^IC)E 



(4) 



The IC losses are given by ^ic = (2l?>) c2B^-^^ and the qu antity 
C2 = 2e^/3mlc^ is a constant defined by (iPacholczyld 1 97Ql) . Two 
main expressions for ^syn exist. 

The Jaff'e-Perola (JP) model takes into account that the pitch 
angles of the synchrotron emitting electrons are continuously 
isotropized on a tim escale shorter than the radiative timescale 
(lJalfe & Perolalll973h 

^syn = (2/3)C2B^ . (5) 

For the Kardashev-Pacholczyk (KP) model, the pitch angle of 
the electrons remains in its original orientation with respect to 
the magnetic field 

^,y^ = C2B^sm\e). (6) 

The JP model is more realistic from a physical point of 
view, as an anisotropic pitch angle distribution will become more 
isotropic due to changes in the magnetic field strength between 
difi'er ent regions and scattering by self-induced Alfve n waves 
re.g. JCarilliet"aDll991l:lTVihbig The JP 

and the KP models assume a single impulsive injection that pro- 
duces a power-law distribution of relativistic electrons. At ^ = 0, 
N(E) = NqE'^ and at some later time (0 



N(E,e,t) = 







if ^ < l/^t. 



(7) 



With this equation for N(E, 0, t) and Eq.[3]the synchrotron spec- 
trum can be computed. The result is that the energy losses (both 
for the JP and KP models) cause the radio spectrum to steepen 
(see Fig.[TT]) above a break frequency (vbrk), with 



B 



Vbrk 



(8) 



JP- spectra have an exponential high-frequency cutoff' above Vbrk- 
KP-spectra steepen to a power-law with a slope (4/3)cKinj - 1 
above Vbrk, because there are always some high energy electrons 
with low pitch angles to radiate at high frequencies in the KP 
case. 

In addition to the above mentioned JP and KP models, where 
all particles are injected in one single burst a nd then simply age , 
there is the continuous injection (CI) model (iPacholcz^l 1 97Ql) . 
i.e., a fresh supply of particles is i njected continuously and 
the so-called KGJP or KGKP models (iKomissarov & Gubanovl 
1994), where particles are injected for a fixed period of time 
after which the supply of newly injected electrons is switched 
off'. This results in different expressi ons for N{E, 6, t), see eq. 6 
from IKomissarov & Gubanovl (11994 . CI and KGJP/KP spectra 
are basically summations of individual JP and KP spectra with 
diff'erent amounts of spectral ageing. Examples of resulting spec- 
tral shapes used in this work are shown in Fig.fTT] 



5.2. Effect of resolution on surface brightness 

If radiative losses are the dominant eff'ect then all spectra for 
diff'erent parts of the relic should line up at some fiducial low 
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Fig. 8. Left: Fitted spectral index map between 2272 and 147 MHz. The spectral index map was made by fitting a second order 
polynomial in log (S) - log (y) for each pixel in the maps at 2272, 1714, 1382, 1221, 610, 325, 241 and 147 MHz. Contour levels 
are drawn at [1, 2, 4, 8, . . .] x 0.175 mJy beam"^ and are from the L-band image in Fig. [3] The spectral index map has a resolution 
of 20'' X 18'' and pixels below 1.5(Trms were blanked. Right: Fitted spectral index (2272 and 147 MHz) for each pixel in the maps 
at 2272, 1714, 1382, 1221, 610, 325, 241 and 147 MHz. The individual maps were convolved with Gaussians of 35" FWHM and 
pixels with a spectral index error > 0.7 were blanked, contours are drawn as in the left panel. 




Fig. 10. Integrated radio spectra for sources C (left), D (middle), and E (right). The lines are fitted straight power-laws with indices 
given in Tabled 



energy/frequency where radiation losses are not important. The 
spectra displayed in Fig. [13] should therefore all line up at 
low enough frequencies. Part of such a trend is indeed visible. 
However, the spectra with a^Qf = -0.65, -0.75 have a lower sur- 
face brightnesses than we naively expect. We attribute this to the 
fact that the pixels, from which these spectra are determined, are 
all located on the sharp outer edge of the relic, where the emis- 
sion drops abruptly to zero within a single beam element. Since 
the resolution of the spectral index map is only about 20" and 
there are about 4 pixels per beam, the regions with the flattest 
spectral index are all located on pixels where the surface bright- 
ness drops at the front of the relic due to the finite beam size. This 
can be seen by looking at the black contours in Fig. [HI i.e., the 
surface brightness does not peak at the region with the flattest 
spectral index. Comparing these contours to the higher resolu- 
tion images in Figs. [H and [6l we find however that the outer rim 
of the relic is indeed the brightest part of the relic. The lower sur- 
face brightness for the Ofref = -0.65, -0.75 regions is thus caused 
by the low spatial resolution of the spectral index map. This is 



also confirmed by the high-resolution 325-610 MHz spectral in- 
dex map, see Fig.[71 

It is interesting to note that the spectra for the regions with a 
flattest spectral index are all relatively straight and show little de- 
viations from a pure power-laws (Fig. [13]). Going to the regions 
with a steeper a^Qf, the spectral shapes change and display signs 
of curvature at the higher frequencies. These changes happens 
gradually from the regions with a flat o^ref to a steeper a^ref . 



5.3. Color-color diagrams 

To better investigate the spectral shape we use so called three 
frequency "color-color" diagrams for which the spectral indices 
for diff' erent region of the sou rce can all be put in a single di- 
agram ( Katz- Stone et a D 11993). On the X-axis we plot the low- 
frequency spectral index and on the y-axis the high-frequency 
one. Color-color diagrams are particularly useful to discriminate 
between theoretical synchrotron spectral models, such as the JP, 
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Fig. 11. Overview of the spectral models used in this work (left) and corresponding radio color-color corves (right) to highlight 
differences in the spectral shapes. All spectral models have Ofinj = 0.6. Solid grey line and filled grey circle (in the color-color 
digram): power-law spectrum without any ageing, solid black line: KP, dotted line: CI, dash-dotted line: KGJP, dashed line: JP. Note 
that a power-law spectrum will translate in a single point on a color-color diagram. The o^iow-frequency = ^^high-frequency line divides 
concave and convex spectral shapes between the frequencies where the spectral index is computed. 



KP, CI, KGJP, KGKP models described before, see Fig. [IH for 
several examples. 

Another advantage of color-color diagrams is that the shapes 
traced out in the diagrams are conserved for changes in the mag- 
netic field, adiabatic expansion or compression, and the radiation 
losses, for standard spectral models. We observe a different por- 
tion of the spectrum for each different frequency or source phys- 
ical condition (e.g., local B field). We can thus use the observed 
curve to constrain the different models (JP, CI, KP, etc.) and in- 
jection spectral indices. Mixing of emission, for example from 
regions with different B strengths or radiation losses, due the fi- 
nite resolution of the observations or projection effects, will lead 
to different curves in the color-color diagrams. If a global spec- 
tral shape exists, it also allows for a better effective frequency 
sampling and provides sampling of a larger range of electrons 
energies. 

5.3.1 . Color-color diagrams for B1 and B2+B3 

The color-color diagrams for Bl and B2-hB3 are shown in 
Fig. [141 (left and right panels, respectively). The points in the 
color-color diagrams seem to trace out single curves. This sug- 
gests the existence of a global spectral shape for these regions. 
The curves are similar for the Bl and B2-I-B3 regions, but they 
slightly differ for a^Qf < -1.4. 

Tracing the curves back to the 0^241 = ^1382 P^wer- 
law) line gives the injection spectral index. In both cases we 
find Qfinj is about -0.6 to -0.7 for Bl and B2-hB3, in agreement 
with the spectral index maps. The color-color diagrams high- 
light the trend of increasing spectral curvature, i.e., distance from 
the power-law line, with decreasing a^cf (Fig. [13]). The spectra 
for Qfref = -0.65, -0.75 are power-laws. We also indicated vari- 
ous spectral models in the diagrams. We compare the observed 
curves in the color-color diagrams to the some of the standard 
spectral models. 

The observed curve is clearly different from the continu- 
ous injection s pectral model. Th e CI model only steepens to 
a = Qfinj - 0.5 (lPacholczvklll970l) . The KP model (with ofinj = 



-0.6) roughly follows the data for flat o^ref, but deviates for 
Q^ref ^ -1. In the end the KP curve bends back and returns to 
the (^24^ = chilli hne. We do not see evidence for a turn back 

to the 0^241 = ^1382 ^^^^ although we do not sample 

this part of the diagram well because of insufficient SNR at high- 
frequencies. As mentioned previously, a KP model is rather un- 
physical. However, it sometimes provides a good fit to observed 
radio spectra because it is similar to the shapes caused by losses 
in a randomly varying magnetic field ( Tribble 1993) . 

The JP model (plotted for ainj = -0.6 and -0.7) follows the 
KP model for flat c^ref . Instead of turning back to the 0^24 ^ = c^HH 
line, the spectral curvature keeps increasing. The data also show 
this trend of increasing curvature, but less quickly than the JP 
curve. A KGJP model with o^inj = -0.7 matches the data for 
B2-I-B3 very well. The KGJP curves in Fig. [141 are for particles 
injected continuously for about 0.6 x 10^ yr and B = 9 yuGauss 
(see Sect. 17.3b . The KGJP model also provides a better match 
for the Bl region compared to the KP and JP curves. Although, 
in the regions with the steepest a^ref values the KGJP model still 
overestimates the amount of curvature. A KGJP model can be 
thought of as an integration of JP spectra with a range of spec- 
tral ages, up to the oldest population of electrons. At first, the 
KGJP model follows the CI model until at some point in time 
the supply of newly injected electrons is shut off. Our injection 
time of 0.6x 10^ yr leads to a distribution of spectral ages in each 
CKref region that is 0.6 x 10^ yr wide at maximum, i.e., it repre- 
sents the amount of mixing encountered in each a^Qf region. It 
should not be interpreted as the total injection time for the relic. 
This value thus depends on the resolution of the spectral maps 
and unavoidable projection effects. 

5.3.2. The effect of resolution and mixing 

One important effect is mixing of emission within the beam. 
Each beam samples superpositions of regions with potentially 
physical different conditions. To investigate the effect of reso- 
lution in the color-color diagrams, we decreased the resolution 
of the spectral index maps by convolving them with Gaussians 
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of 60'' FWHM. The data for the lower resolution maps are also 
displayed in Fig. [141 The effect of the lower resolution is that the 
curve ends up closer to the 0^24^ = c^HH hne. This is expected be- 
cause the total integrated spectrum of the relic has a power-law 
shape (Fig. [9]). 

We also investigated for B 1 , the effect of a possible under- 
lying flux component from the radio halo C. For this, we deter- 
mined the average spectrum of the radio halo and subtracted this 
flux contribution at each frequency. This however did not sub- 
stantially change the resulting color-color diagram for Bl be- 
cause the surface brightness of the relic is much higher than that 
of the radio halo. 

Our conclusion is that mixing of emission (from regions with 
different 5-fields or different electron energy distributions, spec- 
tral ages, etc.) pushes the spectra closer to power-law shapes. 
The spectral curvature we find can thus be regarded as a lower 
limit on the actual curvature and it is therefore important that 
when searching for spectral curvature one retains sufficient spa- 
tial resolution (enough to properly resolve the spectral varia- 
tions). 

5.4. Global spectrum 

Since the data points trace out a well defined curve in the color- 
color diagrams, we have attempted to map the flux measure- 
ments onto a single spec tral sh ape using t he ''shift-t echnique" 
described by Katz- Stone et al.l ([1993); Rud nick & Katz- Stone 
(1996); Rudnickl (1200 1). This allows for a much better sampling 
in frequency, and maps out the spectrum over a larger range of 
electron energies. The idea is that the individual radio spectra for 
different a^cf regions each trace some part of the "global spec- 
trum" (or electron energy distribution) of the source, depending 
on the energy losses, and magnetic fields in these regions. The 
shifts (in log (/)-log (v) space) remove the effects of possible dif- 
ferent local magnetic field strengths, electron densities, radiative 
energy losses, and adiabatic energy gains/losses. 

We shifted the spectra, displayed Fig. [131 in log (/)-log (v) 
space and tried to line them all up. This went remarkably well, 
indicating we indeed have a single electron energy distribution 
(or global spectral shape) that is consistent with the spectra for 
each individual region. The two resulting spectra (for Bl and 
B2-hB3) are displayed in Fig. \T5\ As can be seen it allows for 
an almost continuous sampling of the radio relic spectrum over 
about 4 orders of magnitude in frequency (two orders of magni- 
tude in energy, see Eq.Q. 

Similar to the color-color diagrams we find that the low- 
frequency part of the global spectrum has a power-law shape, 
with a = -0.6 to -0.7. At higher frequencies the spectrum steep- 
ens. We also compare the global spectra with some of the stan- 
dard models in Fig. [151 The results are the same as for the color- 
color diagrams, i.e., the best match we find for a KGJP model 
and for Bl the KGJP spectral curvature is somewhat too high 
for VeflF > 20 GHz. 

5.4.1. Shift diagrams 

Interestingly, the shifts made to align up all the individual ra- 
dio spectra to create Fig. [151 provide inf ormation about changes 
in the underlying physical parameters (iKatz-Stone & Rudnickl 
11994^. The shifts made in the frequency direction (log(y)), to 
line up the spectra, are related to y^B. The shifts in log (/) are 
related to NjB, where A^t is the total number of relativistic elec- 
trons in the volume determined by the beam size. 
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Fig. 12. Schematic representation of two radio spectra with dif- 
ferent spectral ages (break frequencies). To align the red spec- 
trum to the black one in log-log space, a shift in the frequency 
direction is required to match up the break frequencies. The red 
spectrum then falls below the injection power-law spectrum and 
needs to be shifted upwards. The slope A log (I) /A log (v) equals 
the injection spectral index in this case (see the main text). 

By plotting the shifts which we needed to align up the in- 
dividual spectra (log (v) shifts against log (/) shifts), we can in- 
vestigate the contributions of the above mentioned quantities, 
i.e., which quantities cause the spectra to differ from region to 
region. The shift diagrams are shown in Fig. [161 For example, 
a slope (^ log (7)/^ log (y)) of indicates mainly variations in 
the magnetic field strength, a slope of infinity indicates density 
variations, and a slope of zero energy variations. If every line of 
sight samples the same physical conditions (i.e., the same mag- 
netic field, path length through the source, number of relativistic 
particles at some fiducial low energy), then ageing alone would 
give a slope equal to the injection spectral index, see Fig. [121 
This can be seen as follows: one needs to shift the red spectrum 
in Fig. [TJlin log(y) space to match up the break frequency of 
the black spectrum, but then the red spectrum will fall below the 
initial power-law injection spectrum, so a shift (up) in log (/) is 
also required. For the shift in log (v) 

A log (y) = log (ybrk,2) - log (ybrk,i ) • (9) 
For the shift in log / we have 

A log (/) = log (/0Vbrk,2) - log (^OVbru) = ^^^J ^^""^ (^brk,2) " log (Vbrk,!)] (10) 

and thus the slope in the shift diagram is Alog (/)/Alog (y) = 

We find slopes of -0.83 and -0.67 in the shifts diagrams for 
Bl and B2-hB3, respectively. We did not include the first three 
(Bl) and two (B2-hB3) points for fitting the slope. This because 
the log (/)- shifts for these points are affected by the reduced sur- 
face brightness at the outer edge of the relic, see Sect. 15.21 The 
slopes are close to the injection spectral indices of -0.6 to -0.7 
we found from the color-color diagrams. The slope for Bl de- 
viates a little more, but overall the shift diagrams indicate that 
spectral ageing is likely the dominant factor in explaining the 
different spectral shapes from region to region. Apparently, the 
magnetic field and total number of number of relativistic elec- 
trons remain more or less constant for the different Ofref regions. 
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Fig. 13. Normalized radio spectra for relic Bl (top) and B2+B3 (bottom). Each spectrum corresponds to a region from the spectral 
index map where a^ref - 0-05 < a^ref < Q^ref + 0.05, with a^Qf = -0.65, -0.75, -1.85 for Bl, and a^^f = -0.75, -0.85, . . . , -1.65 
for B2+B3. The spectra are plotted in log v - log (v/) space to emphasize the differences between the spectra. First order polynomial 
fits, done in log (/) - log (v) space, are shown by the solid lines to emphasize possible deviations from power-law radio spectra. 
Colors correspond to the different a^ref and go from black, blue, purple, red, to orange. 
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Fig. 14. Color-color diagrams for Bl (left) and B2+B3 (right). Black data points are from the spectra displayed in Fig.[T3l The solid 
black (slightly curved) line is a JP model with ofinj = -0.6, the dashed black line a JP model with ofinj = -0.7, the solid red line a 
KP model with ofinj = -0.6, the dashed red line a CI model with ofinj = -0.6, and the blue line a KGJP model with ofinj = -0.7. 
The grey data points are for the same data as the black points, except the maps were smoothed with a 60'' FWHM Gaussian. The 
^241 ~ ^1382 ^^^^ ^^^^ shown for reference. 




Fig. 15. "Global" radio spectra for Bl (left) and B2+B3 (right). The spectrum for each a^Qf spectral region from Fig. [13] has been 
shifted in log (/) - log (v) space to create the "global" spectral shape. Fourth order polynomial fits through these data points are 
shown by solid thin black lines. JP (red), KP (black) and KGJP (blue) models are shown by the dashed lines. The color coding, 
based on the a^Qf regions, is the same as in Fig.[T3l 



6. RM-synthesis & polarization maps 

6.1. Polarization maps 

The distribution of the polarization E- vectors at 4.9 and 
1.382 GHz are shown in Fig.[T7J In the 4.9 GHz map, we find 
the relic to be polarized over its entire length (in the region with 
sufficient SNR). The E- vectors are mainly perpendicular to the 
relic's orientation, except at the eastern end of Bl (where the 
relic bends and is connected to the linear eastward extension B2). 



The polarization fraction across the relic varies between 10 and 
60%, where it can be measured. At the front of Bl, the polariza- 
tion fraction is mostly between 15 and 30%, while at the front of 
B2 the fraction is 50% or higher. In the 4.9 GHz map we cannot 
determine the polarization fraction for B3 because the attenua- 
tion of the primary beam . It should be noted that the spatially 
averaged polarization fractions are lower. 

At 1382 MHz the polarization fraction drops significantly for 
Bl and B2, compared to at 4.9 GHz. For region B3 we measure 
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Fig. 16. Shifts in log (/) - log (v) space, for Bl (left) and B2+B3 (right), needed to line up the radio spectra displayed in Fig. [15] A 
line is fitted through the data points, excluding the first three (Bl) or two (B2+B3) points. Color coding for the a^^f regions is the 
same as in Fig.[T3l 



a polarization fraction as high as 40%, while for the brightest 
part of B2 it is about 5%, and for Bl it drops below 1%. The 
average polarization fractions and depolarization properties for 
relic B are described in Sect. 16.31 



6.2. RM-synthesis 

Faraday rotation changes the intrinsic polarization angle O^o) de- 
pending on the wavelength (A) or frequency of the radiation. 
The Faraday depth (0) is rela ted to the properties of the plasma 
causi ng the Faraday rotation (iBurnll 19661: iBrentjens & de Bruynl 
[2005h 

-►telescope 



0(r) = 0.81 



HqB • dr [rad m ] 



(11) 



where dr is an infinitesimal path length in parsec along the line 
of sight, B the magnetic field in yuGauss and Hq the electron den- 
sity in cm"^. The sign is defined to be positive for a magnetic 
field pointing towards the observer. The rotation measure (RM) 
is defined as 

If there is only one source along the line of sight (without 
internal Faraday rotation), the Faraday depth is equal the rotation 
measure (RM) at all wavelengths. In other words, all polarized 
emission is observed at a single Faraday depth 0. The observed 
polarization angle (x) is then 



(13) 



In more complicated situations th is relation is not valid (e.g., 
|ValledLl980; Sokoloff'et aUll998|). By expressing the polariza- 
tion as a co mplex vector P = pe^^^ = Q + iU, with p the intrinsic 
polarization lBurnI (Il966h showed that 



-f 



F{(P)e^''^^' d(P 



(14) 



where F((f)) is the Faraday dispersion function, i.e., the complex 
polarized surface brightness per unit Faraday depth. Eq.[T4lis in- 



vertible if the intrinsic polarization angle is constant as func- 
tion of A, then F((p) can be found by measuring P(A^) 



F(cp) 



X+oo 
P(A^)e-^''^'^ dA^ . (15) 
oo 



With modern correlator backends, such as the one at the WSRT 
or EVLA, P(A^) can measured over a large number of frequency 
channel s. If (l)SA^ <^ 1 , with SA^ the ch annel width in wavelength 
squared. ferentj ens & de Bruynl (l2005b showed that Eq.[T5]can be 
approximated as 



F(cp) 



N 

KY^^jPje 



with Wj being some weights and K a normalization factor 

1 



K = 



(16) 



(17) 



As explained by Bren tjens & de Bruynl adding the term /1q in 
the exponent of Eq. [16] results in a better sidelobe pattern for the 
rotation measure spread function (RMSF), with 



and the RMSF given by 
RMSF(0) 



wjAj 



ZN 



(18) 



(19) 



Because P is not measured for every possible Aj, the true 
Faraday depth function is related to Eq.[T6]by a convolution 



F(4,) = F(4,) * RMSF((/i) . 



(20) 



Depending on the sidelobe structure of the RMSF (which can 
be adjusted by choosing the weights wj) a brighter compo- 
nent can contaminate the response of fainter components in 
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Fig. 17. Left: WSRT 4.9 GHz polarization E- vector map. Total polarized intensity is shown as grayscale image. Vectors depict the 
polarization E- vectors, the ir length represents the polarization fraction. The length of the E- vectors are corrected for Ricean bias 
(IWardle & KronbeTglll974^) . A reference vector for a polarization fraction of 100% is shown in the bottom left corner. No vectors 
were drawn for pixels with a SNR < 3 in the total polarized intensity image. Contour levels are drawn at [1, 4, 16, 64, . . .] x Scrims 
and are from the Stokes I image. The beam size is 7.0'' x 4.7'' and indicated in the bottom left corner of the image. Right: WSRT 
1382 MHz polarization E- vector map made with natural weighting. Total polarized intensity is shown as grayscale image. No 
vectors were drawn for pixels with a SNR < 4 in the total polarized intensity image. Contour levels are drawn at [1, 4, 16, 64, . . .] x 
0.18 mJy beam"^. and are from the Stokes I image. 



F((p) and therefore deconvolution might be necessary. The one- 
dimensional deconvolution algo rithm which has been used is 
a simple extension of HogbomI (11974)) RM- CLEAN alg orithm to 
the complex doma in as described by lBrentjensI (l2007l) . see also 
iHeald et al.1 (l2009h . It works as follows 

1. For a spatial pixel the maximum of \F((p)\ is found and the 
location of the peak is stored (0peak). 

2. If this maximum value is higher than a user defined cut- 
off', the RMSF is shifted to the location of the peak (^peak) 
and a scaled version (i.e., using a gain of 0.1) is subtracted 
from F((p). Since the RMSF is complex, a multiplication 
with a phase factor is needed which depends on the phase 

of ^(0peak). 

3. The complex scaling factor, which was used to shift and sub- 
tract the RMSF, is stored as a clean component. 

4. Steps 2 - 4 are repeated until a user defined threshold is 
reached or the number of iterations has reached a predefined 
maximum. 

5. Optionally, the clean components are restored with 
Gaussians with a FWHM matched to that of the RMSF. 
The Gaussians are again shifted to their respective locations 
(0peak), and since this happens in the complex domain, mul- 
tiplied by a phase factor (i.e., the complex part of the clean 
component). The restored Gaussians are added to the resid- 
ual from step 4. If there are multiple RM components within 
the beam it is however not a good idea to reconvolve the 
spectra with the RMSF if one wants to investigate the de- 
polarization properties. This is related to the fact that the 
clean components are co mplex numbe rs and can interfere 
with each other (.Earns worth et al.ll201 ih . 



6. The algorithm then continues with the next spatial pixel. 

The above described RM-syn thesis technique works better 
than linear fitting techniques (e.g.,'Ruzm aikin & Sokoloffll 19791: 
FDolag et al. 2005; Vogt et al. 2005), which break down in low 
signal to noise regimes and do not apply to more complex situa- 
tions when Eq.[T5]is not valid. 

The RM-synthesis technique (i.e., Eq.O and RM-CLEAN al- 
gorithm were implemented in IDL. 



6.2.1 . Application to the L-band WSRT data 

We performed RM-synthesis on the WSRT 25, 21 and 18 cm 
data using cubes of Stokes Q and U images with a resolution 
of 40" X 28", see also Sect. [2]2l These data give a sensitiv- 
ity to polarized emission up to a maximum Faraday depth of 
I0maxl ~ ^ 8.8 X 10^ rad m'^ (Brentie ns & de Bruvnl 

2005). The three bands provide a resolution of ^0 2 V3/A/l^ = 
88 rad m'^, with AA^ = /l^^^ - /l^.^ the total bandwidth in 
wavelength squared. The maximum scale in space to which 
the sensitivity has dropped to 50% is approximately 7i/A^.^ = 
111 rad m"^. The RMSF is shown in Fig [181 Since the first 
sidelobe of the RMSF is about 65% of the main lobe, we 
used the RM-CLEAN algorithm, cleaning down to a threshold 
of 0.12 mJy beam"^ RMSF"^. The distribution of the peak of 
|F(0)|, the so-called "rotation measure map" (or Faraday depth 
map), is shown in Fig.[T9l 

We extracted the location of the peak in the Faraday depth 
(0peak) for all compact sources within 0.5° radius of the clus- 
ter. In total for 17 sources we find an outlier-resistant mean of 
-hi 1.7 ± 3.0 rad m"^, clipping values more than two standard 
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Fig. 18. RMSF from the 25, 21, 18 cm bands WSRT data. The 
solid line is the overall amplitude of the complex RMSF. Dotted 
line depicts the Real part of the RMSF and the dashed line the 
Imaginary part. The FWHM of the RMSF is 88 rad m"^. 




Right Ascension 



Fig. 19. Faraday depth value of the peak in F((p). Contours 
are from the WSRT 1382 MHz image and drawn at levels of 
[1, 3, 9, 27, . . .] X 0.225 mJy beam"!. 



deviations away from the median, with the uncertainty the stan- 
dard deviation of the mean of the 0peak distribution. This is in 
agreement with the value of 1 1 rad m"^ from iTay lor et al.l (l2009l) 
derived from NVSS data measured over scales of 8°.The[F(0)| 
spectra for seven compact sources are displayed in Fig. [2Ql (left 
panel). 

For the relic, 0peak values vary spatially from about -53 to 
+55 rad m"^. For B3, 0peak is about 10 rad m"^, similar to the 
compact sources in the field. For B2, we find a strong gradient 
from -10 to +55 rad m"^ from east to west. For Bl, 0peak varies 
between -10 and -53 rad m"^. The "clean component" spectra 
are clearly resolved, see Fig. [20l (right panel). 

Interestingly, the 0peak values for B 1 and B2 deviate from the 
average galactic foreground of +11.7 + 3.0 rad m"^, indicating 
that some of the Faraday rotation is probably caused by the ICM. 
The 0peak value for B3 indicates it is located in the cluster out- 
skirts and the line of sight towards it does not pass deep into the 
ICM. B2 and Bl seem to be located progressively deeper into 
the ICM of the cluster (or more behind it) given the larger ^peak 
deviation from the galactic foreground. Although the number of 
counts is low in the ROSAT image (Fig. (3]), it is indeed expected 

^telescope 

that Hq dr increases from B3, B2, to B 1 . 

Jsource ^ 

6.3. The depolarization properties of the radio relic 

Our observations show a significant decrease in the polariza- 
tion fraction towards lower frequencies. There are several fac- 
tors that can cause such an eff'ect. We rule out the eff'ect of band- 
width depolarization, as the frequency resolution of the observa- 
tions is sufficient to detected polarized emission with a Faraday 
depth > 10^ rad m ^ using RM-synthesis. A second possibility 
is beam depolarization due to multiple RM components within 
the beam (spatially varying). These diff'erent RM components 
can be caused by variations of the magnetic field and/or electron 
density. The spatially averaged polarization fraction over B 1 and 
B2 is about 13% at 4.9 GHz, 1.6% at 1714 MHz cm and < 1% 
at 1382 and 1221 MHz. We excluded B3 from the analysis be- 



cause the SNR is too low at 4.9 GHz due to the primary beam 
attenuation. 

Two types of depolarization, internal and external, have been 
identified (Burn 1966; Tribble 1991). Internal depolarization 
takes place within the radio source itself, while external depolar- 
ization takes place in a medium between the source and the ob- 
server, e.g., the ICM or magnetized plasma from our own galaxy. 
A simple case of internal depolarization is for a uniform mag- 
netic field and electron density within the source. In this slab 
model (lBurnl il966). the fractional polarization as function of 
/\} is given by a sinc-f unction. From this, the extent in Faraday 
depth of the source can be estimated. Observations indicate that 
the uniform slab model does not fit the data for the radio relics 
in A2256 and the Coma cluster ( Brentjens 2008; Pizzo 201(J). 

For Bl and B2 there is a clear trend of depolarization to- 
wards higher A^. The only way to get depolarization as a function 
of wavelength is to have multiple RM components within the 
beam, either across it or along it - external -, or mixed with the 
synchrotron emitting plasma - internal. Therefore, the Faraday 
spectrum must be resolved wherever there is depolarization. In 
practice, because of the finite width of the RMSF, it may not 
always be possible to see these multiple RM components, even 
though they are depolarizing. In the case of B 1 and B2, the spec- 
tra are indeed clearly resolved, consistent with the high amount 
of depolarization found for these regions. It cannot be directly 
determined whether this is due to beam depolarization (spatially 
varying RM components) or due to multiple RM components 
along the line of sight. In the first case, the depolarization should 
vanish when going to high enough spatial resolution. The fact 
that the RM varies spatially from beam to beam across B 1 and 
B2 makes it very likely that at least a significant part of the de- 
polarization is indeed caused beam depolarization. 

The trend from east to west over the relic of (i) a larger 0peak 
deviation from the galactic foreground, and (ii) an increasing de- 
polarization towards higher suggest that most of the depolar- 
ization is probably due to the ICM of the cluster itself and not 
due to internal depolarization. 
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Fig. 20. Left: Gaussian restored \F((f))\ spectra of compact sources. B3 0559+422A black solid line, NVSS J060406+422125 blue 
dashed line, NVSS J0604 16+4221 10 (south component) blue solid line, NVSS J0604 16+4221 10 (north component) dashed-dotted 
blue line, NVSS J06023 1+42281 1 (west component) red dashed line, NVSS J06023 1+42281 1 (east component) red solid line, and 
NVSS J060206 +422045 black dashed line. The high "residuals" for B3 05 5 9 +422 A are an artifact due to remaining calibration 
errors. Right: \F((p)\ spectra, consisting of clean components, extracted at three different positions from relic B. The positions are 
indicated by the circles drawn in Fig.[T9j Red line: position B3, blue line: position B2, and black line: position Bl (see also Fig. [3]). 



7. Discussion 

The complex diffuse radio emission and extended X-ray emis- 
sion give strong support for the fact that the cluster is cur- 
rently undergoing a merger event. There are however a few puz- 
zling aspects to the diffuse emission in this cluster. The lin- 
ear extension of relic B to the east is quite peculiar. Usually 
shock surfaces are curved a nd not very linear over distances 
of ni o re than a Mpc (e . g.. iHoefteta l.' '2008'; 'BattagU a et"aD 
120091: IVazza et al.l 120091: ISkillman et al. .2011.) . The ongoing 
merger event for IRXS J0603.3+4214 is probably more com- 
plex than the "simple" binary merger events that are tho ught to 
give rise to symmetric double radio relics (e.g., Roettiger et al.l 
[199^;[ van Weeren et al. 201 la). Deep X-ray observations will be 
needed to investigate whether the relic is indeed associated with 
such a linear shock front. 

Relic E could also trace a shock wave with DSA, given its 
large extent and relatively flat integrated spectral index. The na- 
ture of relic D remains unclear, its small size and morphology 
would suggest a radio phoenix scenario, but the integrated spec- 
tral index of -1.10 + 0.05 is not consistent with this interpreta- 
tion. Higher-resolution observations (~ V) would be helpful to 
better delineate the morphology. 

An interesting point is the relation between relic B and the 
radio halo C. The radio halo is located directly south of Bl, the 
brightest part of relic B. The halo and relic are connected by 
a region with a < -2, but then the spectral index flattens to 
a ~ -1.2 and increases a bit towards the center of the halo. 
We speculate that relativistic particles, previously accelerated at 
the shock front, are re- accelerated due to merger induced tur- 
bulence and then form the radio halo. Another interesting fea- 
ture of the halo is the increased surface brightness at its southern 
end. Why does the surface brightness increase here? For a binary 
cluster merger event one expects two shock waves to fo rm travel- 
ing diametrically outwards (e.g. /Roettiger et aO 19991) . If relic B 
traces the main shock wave one could also expect a second shock 
waves at the south side of the cluster center. This could be the 
region of increased surface brightness at the southern end of C. 



Although, this remains somewhat speculative, also because the 
merger event seems to be complex given the peculiar linear ex- 
tension of relic B. X-ray observations could reveal whether the 
south part of radio halo C is associated with a shocked region 
with a higher ICM temperature. 

There is also a possibility that source C is not a radio halo, 
but another relic seen close to face-on, because the surface 
brightness does not clearly peak at the location of the ROSAT X- 
ray peak. However, the radio power of source C is in agreement 
with the Lx-PiAGUz correlation. Also, the brightness distribution 
over C is regular and quite smooth, which argues against a relic 
interpretation. Relics are often compos ed of smaller filamentary 
structures (e.g. JClarke & EnBlm 2006). We note that source C is 
somewhat similar to the diffuse emission s een be tween the relics 
in CIZA J2242.8+5301 dvan Weeren et al.ll2010l) . 



7.1. Relic spectra 

The analysis of the radio spectra for radio relic B are consis- 
tent with the existence a two global spectral shapes, one that de- 
scribes the spectrum for the east part and one for the west part. 
This because it is possible to line up all radio spectra onto a sin- 
gle shape with shifts along log (/) and log (v). 

The spectral shapes of the west (Bl) and east (B2+B3) 
regions only differ in regions in the spectral index map with 
a < -lA (e.g., see Fig.O. At the front part of the relic we find 
power-law radio spectra, consistent with an injection spectral in- 
dex of -0.6 to -0.7. The injection spectral index is also con- 
sistent with the total integrated spectral index of -1.1, i.e., the 
integrated spectral index sho uld be -0.5 units steeper than the 
injection spectral index (e.g.. iMiniatil 120021: iBagchi et al.ll2002h 
for a simple shock model (with the properties of the shock not 
changing over time). We also find that the shape of the radio 
spectra are mostly governed by the effects of spectral ageing, by 
looking at the shifts made to align up all individual spectra to the 
global spectral shape (see Sect. 15.4.1]) . Changes in the magnetic 
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field, total electron content, or electron energy (except energy 
changes due to ageing) do not seem to play a major role. 

The observed radio spectra are not consistent with a JP 
model (which consists of a single burst of injection and spectral 
ageing). Instead, we find that the KGJP model, which consist of 
a sum of individual JP spectra with diff'erent amounts of ageing, 
provides a good match to the spectrum of the east part of the 
relic. This can be thought of as mixing of emission, most likely 
because of projection eff'ects. For the west part of the relic (in 
the regions with the steepest spectral index), we find deviations 
from the KGJP model. We speculate that additional mixing of 
emission takes place here, or the spectra are contaminated by the 
emission from radio halo C (although we attempted to subtract 
this emission). 

The KGJP spectra indicate that (1) at the front of the relic we 
see freshly injected radio plasma, giving a power-law spectrum. 
Little mixing of emission takes place here. (2) For the middle 
part of the relic (the region -0.9 < a < -1.1 from the spectral 
index map), we observe a mix of radio aged radio plasma and 
freshly injected radio plasma. For the back part of the relic (the 
region with a < -1.1), we see a mix of aged radio plasma, but 
no freshly injected plasma. 

Mixing of emission thus seems to be an important factor 
determining the spectral shape. This suppresses spectral curva- 
ture and pushes back the radio spectra to power-law shapes. It is 
therefore crucial that this eff'ect is considered. The lack of spec- 
tral curvature does therefore not directly imply spectral ageing 
is not important. 

An injection spectral index of -0.6 in combination with 
DSA, implicates a Mach number between 4.6, quite high. The 
minimum allowed injection spectral index of -0.7, gives a Mach 
number of 3.3. Instead of DSA, another possibility is shock re- 
acceleration of fossil electrons. In practice, t hese radio spectra 
will probably be indistinguishable from DSA (iMarkevitch et al.l 
120051) . The At > 3 we find is interesting. Simulations indi- 
cate that merger shocks with these Mach nu mbers are rare (e.g., 
IVazza etal.ll201ll: lArava-Melo et al.ll2012h . Deep all-sky radio 
surveys are therefore very important to pinpoint these extreme 
cluster merger shocks. 

7.2. Alternative models to explain the relativistic electrons 
from the radio relic 

An alternative explanation for radio relics, based on a sec ondary 
cosmic ray electron model, has been proposed by lKeshetl (|2010). 
This model predicts a spectral index of -1 at the front of the 
relic. We do not observe this, although there is a hint of a steeper 
spectral index at the front of Bl at some places in the 325- 
610 MHz spectral index map. However, this could be an artifact, 
as edge pixels is the spectral index map are not reliable. Also 
these steeper regions are smaller than the beam size. We do not 
find indications of magnetic field changes across the relic from 
the shifts needed to align up the spectra from individual regions. 
On the other hand, in the images we see a morphological connec- 
tion between the radio relic an d halo, these h alo-relic "bridges" 

are included in the model from Keshefl (120101), 

EnBlin & Gopal-Krishna (200l|) and lEnBlin & BriiggenI 
(I2002D proposed that (some) radio relics could be the result of 
adiabatically compressed fossil radio plasma due to a shock 
wave. This could boost the radio luminosity by about two or- 
ders of magnitude. The question is whether the large relic in 
IRXS J0603.3-h4214 can be exp lained in this way. The model of 
lEnBlin & Gopal-Krishnal ([2QQll) consist of five phases: (1) injec- 
tion by a radio galaxy, (2) expansion of the radio cocoon, during strength scales with (1 + ^) 



this phase the cocoon becomes undetectable, (3) the "lurking"- 
phase, when pressure equilibrium is reached and the volume of 
the cocoon remains constant. Due to significant adiabatic losses 
most of the electrons reside at low energies which reduces the 
radiation losses. The losses can be further reduced by a low- 
magnetic field. (4) Adiabatic compression of the cocoon, the ra- 
dio emission is boosted so it becomes visible again and the break 
frequency moves upwards. (5) The radio emission fades away 
due to synchrotron and IC losses. 

In the case of IRXS J0603.3-h4214, the fossil plasma should 
be quite old (> 1 Gyr) as it takes a considerable amount of time 
for a radio galaxy to move over a distance of 2 Mpc and dump 
the radio plasma. The ch ange in the break frequ ency from the 
adiabatic compression is (iGiacintucci et al.|[2008l) 



5Af-l 



(21) 



Adopting a reasonable val ue of v^^^f = 50 MHz for 1 Gyr 
old fossil radio plasma (e.g.. lEnBlin & Gopal-Krishnal200 ih , we 
would need a At > 8 shock to push y^^^^ above 4 GHz. For more 
realistic Mach numbers of Al < 4 we should thus observe steep 
and curved radio spectra. This contradicts the flat spectral in- 
dices of (2 = -0.6 we find and the straight power-law integrated 
radio spectrum. 

Only in the unlikely scenario that the plasma is about 0. 1 Gyr 
old, before it is compressed by the shock, the break frequency 
could be pushed up to a few GHz. We also note that to produce 
a spectral index gradient a high mass load is required, which 
would allow the shock wave to successively compress the ra- 
dio plasma. In the case of a high mass load due to undetectable 
cool gas, a shock forms and the adiabatic gains are restricted 
due to the limited compressio n factor of such a shock wave 
(lEnBlin & Gopal-Krishnall200lh . 



7.3. Magnetic field 

The RM for the west part of the main relic (B 1) deviates from the 
foreground RM, determined from compact sources in the field 
around the cluster. This suggests that the west part of the relic 
is located deeper into the ICM, also consistent with the higher 
amount of depolarization found for that part of the relic. A con- 
stant magnetic field of 1 yuGauss, = 1 x 10"^ cm"^ and a path 
length 500 kpc, will give a Faraday depth of about 40 rad m"^. 
The observed diff'erence of about 60 rad m"^, compared to the 
galactic foreground value, suggests a relatively high magnetic 
field (B > I yuGauss). Also because the magnetic field topol- 
ogy is likely more complex, i.e., the magnetic field could have 
reversals along the line of sight. 

We estimate the magnetic field at the front of relic Bl 
by assuming minimum energy densities in the relics (e.g., 
iGovoni & Feretti 2004 ) and using the same procedure as de- 
scribed in van Weeren eLalJ ([2009). We take k = 100, i.e, 
the ratio of energy in relativistic protons to that in electrons, 
250 kpc for the depth (d) along the line of sight, a spectral in- 
dex of -0.6, and a surface brightness of 0.205 mJy arcsec"^ at 
610 MHz. We use low and high energy cutoff's (ymin and ymax, 
with ymin ^ Tmax) instead of low and hi gh frequency cutoff's 
(IBeck & Kraus3l2Q05l:lBrunetti et al.lll997h . giving the so-caUed 
"revised" equipartition magnetic strength (B^^). For ymin = 100 
this gives B'^^ = 9.2 yuGauss. If we take ymin = 5000 we get 



^eq - /^Gauss. The revised equipartition magnetic field 
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Both the polarization and equipartition results indicate a rel- 
atively high magnetic field strength. These relatively hig h val- 
ues ha v e also been foun d for ot her relics by Finoguenov eLalJ 
(l2QlQh : Ivan Weeren et al.1 (|2Q1Q|) using limits on the IC X-ray 
emission or modeling the relic's brightness profile. 

8. Conclusions 

We presented detailed GMRT and WSRT radio observations of 
a newly discovered cluster of galaxies at z = 0.225. We find 
the cluster to host three radio relics and also a giant ~ 2 Mpc 
radio halo. The radio power follows the Lx-Piaguz correlation 
for giant radio halos. The extended X-ray emission and complex 
diff'use radio emission indicate we are witnessing an ongoing 
cluster merger event. 

The WSRT observations reveal that the front of the relic is 
highly polarized, with a polarization fraction of up to ~ 60% 
at 4.9 GHz. At lower frequencies (in the L-band), the polar- 
ization fraction drops considerably , as has been se en for other 
well studies relics , e.g., in A2256 (Brentjens 2008) and A2255 
(iPizzo et al.l l201ll). The observed depolarization and Faraday 
depth suggest that the west part of the main relic is located 
deeper into the ICM of the cluster. Using equipartition argument 
we find a high magnetic field strength of about 7-9 yuGauss for 
the bright relic. 

For the radio halo we find a spectral index of a = -1.15 + 
0.06. The spectrum for the bright radio relic is a power-law be- 
tween 74 MHz and 4.9 GHz with a = -1.10 ± 0.02. If this relic 
traces a shock, where particles are accelerated by the DS A mech- 
anism, then the injection spectral index of -0.7 to -0.6 would 
imply a Mach number of 3.3 to 4.6. The spectral index steepens 
systematically across the width of the relic from about -0.6 to 

2. For the bright western part of the relic we can trace the 

spectral steepening to a < -2.5. The radio spectrum at the front 
of the relic has a power-law shape with a = -0.6 to -0.7, while 
the amounts of spectral curvature increases gradually towards 
the back of the relic. 

We analyzed the radio spectra making use of color-color and 
shift diagrams ( Katz- Stone et al. 1993; Katz- Stone & Rudnick 
Il994l) and divided the relic into an eastern and western part. 
Both parts were further divided into smaller regions based on 
the spectral index map in Fig. [8] (left panel). We found that the 
individual radio spectra in each region could all be lined up. This 
implies the existence of a single global electron energy distribu- 
tion. The shifts needed to align up the individual spectra, indi- 
cate that spectral ageing is the dominant factor explaining the 
changes in the spectra from one region to another. Changes in 
the magnetic field, total number of electrons, or adiabatic expan- 
sion/compression do not seem to be important. 

In addition, we find evidence for mixing of radio emission 
with diff'erent amounts of spectral ageing within the beam. The 
amount of mixing increases away from the front of the relic. A 
so-called KGJP model (see Sect. l5.1l) , which can be thought of as 
a mix of spectra with diff'erent ages, describes the global spectral 
shape for the east part of the relic. For the west part of the relic, 
in particular the region with the steepest spectral index, we find 
a small deviation from this shape, possibly because the amount 
of mixing increases here. 

The eff'ects of spectral curvature are suppressed by mix- 
ing and this pushes the radio spectra closer towards power-law 
shapes. This means that if a spectral index gradient is seen, but 
no spectral curvature (or only very little), this does not necessar- 
ily imply that the spectral index gradient is not related to spectral 
ageing. 



We conclude that the color-color and shift diagrams provide 
a powerful tool to constrain physical conditions that shape the 
radio spectra. In the future we plan (i) to expand the color-color 
analyses to smaller regions, by obtaining higher resolution maps, 

(ii) increase the frequency range to < 150 MHz or > 3 GHz, and 

(iii) apply the technique to other radio relics and halos. 

For the origin of the relativistic electrons (in the bright relic), 
the results favor a scenario where particles are accelerated or re- 
accelerated in a shock. Interesting is the high Mach number we 
found, simulations indicate that merger shocks with these Mach 
numbers are rare. A crucial test would be the to measure the 
Mach number directly from the jump conditions and compare 
this with Qfinj . Deep X-ray observations of this cluster and sim- 
ulations are also needed to constrain the merger scenario and 
determine the dynamical state of the cluster. 
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Appendix A: Spectral index error maps 
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Fig. A.l. Spectral index error map corresponding to Fig.lH The map is computed on the basis of (Trms values for the individual maps 
and the uncertainties in the absolute flux scale. 
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Fig. A.2. Left: Spectral index error map corresponding to Fig.[8]left panel. Right: Spectral index error map corresponding to Fig. [8] 
right panel. 
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